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ABSTRACT 

A comprehensive analysis of the Vuilleumier refrigerator was con- 
ducted. This analysis includes the effects of non- isothermal gas heat 
addition and rejection, hot and cold regenerator inefficiencies, conduc- 
tion losses, and gas leakage losses. A computer program was written 
which solves the equations resulting from the analysis. The program 
calculates internal pressures, temperatures, and gas flow rates as func- 
tions of refrigerator crank angle, as well as overall refrigerator cooling 
load and power input. Comparisons between the program results and 
available data show good agreement, with a marked improvement over 
the predictions of the ideal model. 
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MATHEMATICAL ANALYSIS OF A 


VUILLEUMIER REFRIGERATOR 


by 

Allan Sherman 
Goddard Space Flight Center 


I. INTRODUCTION 

The Vuilleumier (VM) refrigeration concept was patented in the United States in 1918. Ap- 
parently, because of its low efficiency and cooling load when compared with other low-tempera- 
ture refrigerators, its development was not pursued. However, there has recently been much 
interest in the development of a VM cryogenic refrigerator for spacecraft application, because a 
refrigerator of this type has the potential advantages of long lifetime operation, compactness, and 
low weight. 

Figure 1 shows a possible design for a 
cryogenic VM refrigerator for spacecraft ap- 
plication. The refrigerator consists of a hot 
cylinder to which heat is transferred from a 
heat source, a cold cylinder to which heat is 
transferred from the load(e.g., spacecraft sen- 
sors), and a sump from which heat is rejected 
to the ambient (e.g., spacecraft radiator). The 
gas (e.g., helium) is at a different temperature 
in each of these three sections (Tj^ > T^ > T^ )*, 
and continuously flows through the regenerators 
into each section, due to the motion of the dis- 
placers as the crankshaft rotates. The gas temperature differences among the sections are main- 
tained by the displacers, which are thermal barriers; the regenerators, which heat or cool the 
flowing gases to the appropriate temperatures; and the heat transfer to or from each section. 


HEAT TRANSFERRED FROM 
REFRIGERATION LOAD 

o 



Figure 1. VM Refrigerator 


*For discussion purposes, we assume that the individual cylinder gas temperatures are constant. Actually, as will 
be pointed out later, the cylinder gas temperatures fluctuate during the cycle. 
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With reference to Figure 1, the reasons for the advantages of the VM refrigerator for space- 
craft application are: 

(1) Compactness and Low Weight: By its very design the VM refrigerator is essentially a 
one-component system (i.e., the compressor and expander are in one housing), and, for the low re- 
quired cooling loads of spacecraft sensors, is of relatively low weight. Furthermore, as will be 
discussed, the power required to drive the VM is very low, and thus a large motor is not required. 
However, a heat source (e.g., a radiosotope type) is needed. 

(2) Long Lifetime Operation: The pressure differences among the three sections of the VM 
refrigerator are very small. Hence, piston pressure seals are not required (as contrasted to other 
possible refrigeration concepts such as the Stirling cycle), and, in fact, a feasible seal design might 
be just a clearance between the displacers and cylinder walls. Thus failure by pressure seal wear 
could be eliminated and lower rolling bearing forces could be attained. 

II. BASIC VM OPERATION 

The VM cooler (Figure 7.) is essentially a heat engine driving a refrigerator. Thus the input 
energy to the device is supplied as heat, which minimizes the required drive- motor size, and allows 
the direct use of a long lifetime heat source such as a radioactive isotope. The displace, ~nove the 
gas from one cylinder of the refrigerator to the other, but do not directly compress or expand the 
gas. 


The steady state operation of the VM refrigeration cycle will now be explained. Reference 
will be made to Figure 2 (a schematic of four crank positions encountered during operation, 
and Figure 3 (the pressure-volume diagrams for the cold cylinder, hot cylinder, and total gas 
volume for the VM cycle) . In this explanation we assume that the pressure drop across the re- 
generators is zero, and thus for the design depicted the pressures in the three sections are 
always equal. (In the actual case, the regenerator pressure drops are comparatively small.) 
^Furthermore, we note that the shapes of the pressure-volume diagrams in Figure 3 are close to 
those resulting from the crank type of design shown in Figure 1. One could conceive of a displacer 
movement scheme which would result in more idealized VM performance (i.e., higher thermal 
•sfficiency), but the required mechanism would not be practical for long time spacecraft use. Finally, 
:it should be noted that the following is a simplified view of the operation of the VM engine, in that 
only the predominant processes for a given crank position are pointed out. The actual operation of 
’:he refrigerator is much more complicated, as the subsequent mathematical analysis will show. 

We begin with the crank in the South position (Figure 2). At this point the cold displacer is 
at its maximum displacement position, and the hot displacer is only at its half- maximum position. 
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Figure 2. VM Refrigerator Operation 



Vh 


The mean gas temperature in the VM refrigera- 
tor is relatively low, and consequently the gas 
pressure is low. This is shown in the hot cyl- 
inder and cold cylinder P- V diagrams in Figure 
3. Since this low pressure is the result of the 
West-South expansion process (as will be shown), 
heat is absorbed from the refrigeration load and 
heat source. 


N-W 


TOTAL 


E-S 


V 


As the crank turns to the East position, Figure 3. Pressure-Volume Diagrams 

both the hot and cold cylinder volumes decrease. Cycle 

Part of the cold gas is forced through the cold regenerator, which is at some mean temperature, 
Tcrg * (X: "^cRG where it is heated to nearly before entering the ambient section. 

Similarly, most of the hot cylinder gas is forced through the hot regenerator which is at some mean 
temperature, (T^ < < T„ ), where it is cooled to nearly before entering the ambient 

section. Heat, which is to be used later in the cycle, is thus stored in the hot regenerator. Finally, 
since both the hot cylinder volume and cold cylinder volume have decreased, the mean gas temper- 
ature, and consequently the gas pressure, change very little during this process. Nevertheless, 


*The temperature along the length of the hot and cold regenerator matrices varies from a temperature somewhat less 
than the gas temperature at the hot end to somewhat greater than the gas temperature at the cold end. 


3 


Kome t^as expansion (caused by the small pressure drop that does occur), with the resulting heat 
absorption from the load, does take place. 

As the crank moves from East to North, the hot cylinder volume increases and the cold 
cylinder volume decreast s. The cold gas which is forced through the cold regenerator is hea^^'d 
<0 nearly , wh.'le part of the ambient gas flows through the hot regenerator, is heated (from 
stored energy) to nearly T„ , and enters the hot cylinder. The net effect of the hot cylinder 
volumetric increase and the cold cylinder decrease is an increase in the mean gas temperature 
and gas pressure, rience, this process is one of gas compression, and for the temperatures to 
remain constant, heat must be rejected. Therefore, heat is rejected at the ambient section. 

As the crank turns from North to West, the volumes of both the hot cylinder and cold cylinder 
increase. Part of the ambient gas moves through the cold regenerator, releases heat to it, and 
enters the cold volume at nearly T^, . On the hot side, part of the ambient gas moves through the 
hot regenerator, absorbs heat, and enters the hot volume at nearly T^. . Since both the cold volume 
and hot volume increase, the system pressure does not greatly change. Nevertheless, there is 
some compression with a corresponding heat rejection at the ambient section. 

The crank now turns from West to South, decreasing the hot cylinder volume while increas- 
ing the cold cylinder volume. Part of the hot cylinder gas is forced through the hot regenerator 
where it is cooled to nearly T,^ , while part of the ambient gas is forced through the cold regenera- 
tor where it is cooled to nearly T^,. The mean gas temperature, and consequently the gas pres- 
sure, decreases, resulting in gas expansion. This causes heat absorption at the cold and hot ends. 


III. ISOTHERMAL IDEAL MODEL 

The isothermal ideal mathematical model assumes isothermal expansion and compression 
of the gas in each cylinder, 100% thermally efficient regenerators, zero pressure drop across the 
regenerators, no displacer sliding friction, and no beari;ng friction. From the first law of thermo- 
dynamics for the hot cylinder over one cycle. 


^ d 0 = d U + d V - hj 


d m, 


For a regenerator with 100% thermal efficiency. 




hj d - 0, 


( 1 ) 


( 2 ) 
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and since t>te <;as in the cylinder returns to its initial state after one cycle, 


I 


cl u 


0 . 


( 3 ) 


Thus, over one cycle, 


Similarly, 


Qh 


I 


p V, 




f 


p d Vj.. 


(4) 


(5) 


Now, from the equation of state and the conservation of mass, we have at any time during the cycle 


P^Mo 




^CRQj 


( 6 ) 


The volumes Vj. » can be computed as functions of the crank angle. V„„, V^, T„, T^., T^ , 

Tcr > and T„k ^ 1*6 constants. Hence, Equations (4) and (5) can be expressed in terms of the 
variable % and integrated from = 0 to ^ = 2 77 (one cycle). Explicit expressions for the cooling 
loadQ^ and hot cylinder power requirement Q„ can then be obtained (Reference 1). 

The inadequacy of this model becomes all too obvious when we compare experimental results 
with theoretical predictions. Typically, the above analysis may predict a cooling load of 10 watts 
and a hot cylinder heat input of 50 watts when the actual valties are about 2 watts and 180 watts, 
respectively. These discrepancies prompted investigators to correct the ideal model by estimat- 
ing various losses, such as regenerator inefficiency, for the whole cycle, and then adding these 
losses to the ideal model results. This process of uncoupling the losses from the engine per- 
formance resulted in only limited success, with little confiden e in the predictions because of a 
lack of theoretical justification for such a process. 

The analysis presented in this report was initiated with the goal of providing a program for 
more confidently predicting VM refrigerator performance, and moi e accurately weighing design 
trade-offs. An attempt was made to include os many losses as possible without the uncoupling as- 
sumptions. In order to do this, all processes were expressed in the form of pai Hal differential 
equations, and the resulting rather complicated system of equations was solved via a computer 
program. 
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IV. ANALYSIS 


A. General Configuration 

Figure 4 shows a schematic of a VM refrigerator with the significant varial)les employed for 
the analysis. It should be emphasized that this figure is only a schematic and that the analysis is 
not limited to the particular geometry shown. For example, the regenerators could be located with- 
in the refrigerator without affecting the subsequent equations. Furthermore, the directions of the 
arrows in the figure are arbitrary. 

The heat exchange area for each cylinder is subdivided into the cylinder wall and heat ex- 
changer entrance regions. In this way, flexibility is obtained in matching the actual VM geometry 
.0 that of the mathematical model. 

The walls of the cylinder and heat exchanger are assumed to be at a constant temperature 
(T„„, T^^ , Tj,^) for each cylinder section. This means that we are as suming the cylinder walls 
and end caps to be infinite heat sinks (or sources) relative to the fluctuating gas temperatures. 

This is indeed close to the actual case for most VM refrigerator designs. 



Figure 4. VM Refrigerator Schematic 
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B. Geometric Equations 


A schematic of the geometric iirrangement 
of the VM refrigerator with dimensions labeled 
is shown in Figur*: 5. The crank angle is 
measured from the axis of the cold cylinder and 
thus at =0 the cold displacer is at top dead 
center. 


From simple geometric considerations, 
then, it is shown in Appendix A that 
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1 - cos 
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+ 1 - cos ( 
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Figure 5. Geometric Arrangement of 
the VM Refrigerator 
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Differentiating Equations (7) through (9) with respect to time, we have 


^ _ ''CMAXD 

Vc - 


r ■ ,n\ oj cos 6 

T sin t? — — 


a; sin ^ 


^ '^HMAXD 


and 




d -V) 




( 10 ) 

( 11 ) 

(12) 
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These geometric equations describe the cylinder volumes and time rates of change of the volumes 
as functions of crank position, crank assembly dimensions, and rotational speed, and are used later 
in the analysis. 

C. Cylinder Equations 

Referring to Figure 4 for the hot cylinder, we have from the First Law of Thermodynamics 
(neglecting gas kinetic energy terms), 


fl Q fl n + - hj (! tiij - H nij 


(13) 


where, from the perfect gas law, 


and we take 


Ph - ~ R T„. 


(14) 


’’l '' ”l = 1'2H 


2H 


d t, 


(15) 


(16) 

Note that in Equation (16) we are assuming that the leakage gas around the displacer enters and 
leaves the hot cylinder at T„ . This assumption is exact when the gas leaks from the hot cylinder 
to the ambient cylinder: it is a reasonable simplifying approximation when the gas flows from the 
ambient cylinder to the hot section. 

Now 


^ = [«HA - V + d t (17) 

and 

d U = d (M„ Cy T„) = Cy d T„ + Cy T„ d M„. ( 18 ) 

Combining and simplifying Equations (13) through (18) and transforming the independent variable 
via 

dS = (19) 
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yields 


HT, 


H 1 


H '' 


clM., 


~ *SlA ^^A “ fj t 


M,. f'v 


+ - 


(y - 1) T„ v„ 


Cy 


M.. 


*'*>« 'Si 


( 20 ) 


where 


. 

d t ■ *'*2 H ^ ’'’hL- 


(21) 


From a similar analysis for the cold and ambient sections, we have 


H t 


c 1 


Where 


and 


d M. 


»CA - Tc) »^CA (Ta - Tg) ^c— (y - n Vc y m, 


Me Cy 


'CL 


McC, 


M. 


M, 




cl a 


d M. 


“aa^T^w-V K„,(T„-Te) Ke* (T, - Te) T, — (y - 1) T,, V, 


L 




MaC, 


^<hTa 2 >T,Mee y T, M„, 


M. 




n M + M 

ft t 2C + ""cL’ 


d M. 




( 22 ) 


(23) 


(24) 


(25) 


D. Heat Exchange Sections 


For the heat exchange sections, we assume the total volume and hence the gas storage to be 
negligible. Thus we use the usual steady state heat exchanger equations for a constant wall tem- 
perature. We can then say 


^HOUT = \ I - e 


,w c 


(I"* “ "1'hin) ■*■ 1^1 


‘HIN 


(26) 


where the variables are assigned as shown in Table 1. We also assume that the pressure drops 
across the heat exchange sections are negligible compared to the regenerator gas pressure drops. 
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Table 1 


Variables Used In the Heot Exchanger Equation (26) 


Cold Heaf Exchanger 

Hot Heat Exchanger 


P > 

o 

* CRG 

P < P 

C CRG 

p 

II 

< P 

HRG 


> P 

- ‘ hPG 

1 

h*‘ 

‘^CHR 

h** -- H 

CHE 

h** 

H 

MHE 

h** 

^IHE 


T 

Miin 

T 

* 2C 

T ^ T 

HIN ‘ic 

T 

HIN 

T 

IH 

T 

HIN 

T 

'2H 


T 

HOUT 

T 

MC 

^HOUT ^2C 

T 

HOUT 

' 2H 

T 

*HOUT 

T... 


W 

M,c 

W = M^c 

W 


W 



T - 
w 

T 

T T 

w 'cw 

T 

’ w 

' T,iw 

T 
^ w 

T„. 

1 Ambient hieat Exchanger Number 1 J 

Ambient Heat Excha 

nger Number 2 



P 

CRti 

*^A ^ ^CRG 

Pa 

< P 

HRG 


- ^HRG 


h** = 

*^AHEI 

h*’ ^ »A„E, 

h” 

' ^^AHE2 

h“ 

*^*AHE2 


T 

T 

*oc 

Sun *ai 

T 

*HIN 

- T 

OH 

T 

HIN 


1 

1 

I 

T = 

HOUT 

Ta, 

T T 

Nout oc 

T 

HOUT 

- 

T 

*HOUT 

Ton 

1 

w = 


W 

w 


w 



T 

w 

Taw 


T 


T 



E. Regenerators (Hot or Cold) 

For the purposes of this analysis, we assume that the regenerators are of uniform geometry 
and material along their lengths, and describable by a free-flow area and a wetted perimeter 
P*. In addition, we are assuming that both the gas in each regenerator and the material of the re- 
generators have constant specific heat. 

Assuming me-dirnensional flow, the energy equation (dropping the gas kinetic energj' term) 
for an infinitesriual slab of gJis, perpendicular to the regenerator axis, and of mass dm is 


d 0 
d t 


9 h 


C X 



Now, 


^ P « B m 

— — — + 1 
9x 


(27) 


( 28 ) 
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r 



and 


u 

TF 


>") Cy 


t 


T 

g V 


i (H m) 
~ t 


v\. ''re 


dm ■ / rl X. 


In addition we may write 


and 


c' m A ^ / 


^=h*P‘(T,-T^)dx. 


(29) 

(30) 


(31) 

(32) 


Combining and simplifying Equations (27) through (32), while changing the independent vari- 
able from t to yields 


S _ 1 ^ (7 - 1) h* P* (T, - T^) - m (y - D J (-^J - R (7 - 1) 


B 6 ^ Oy 


d m 

9 x4 




(33) 


Now, an energy balance in the corresponding meta_. regenerator slab, ignoring axial conduc- 
tion, is 


Llk _ 1 h* ri T 'i 


(34) 


from which 




1 r h* 

K/Lro) 


"TR 


(35) 


The convective heat transfer coefficient h in Equations (32) - (34) is assumed to be of the 


form 


h*=DmE, (36) 

where D and E are obtained from experimental data (e.^.. Reference 3). 

In order to solve Equations (32) and (35), and to calculate the corresponding mass flow rates, 
a ’’lumped" type of approximation similar to that used in Reference 2 is employed. For the purpose 
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Figure 6. Illustrated of Lumped Flow 
Resistance Approximation 


of this approximation the gas in each regenerator is considered to be at an average temperature 


T given by 




(x) fl X, 


(37) 


and an average pressure given by 


M 

R T . 

v « 

RG 


(38) 


The flow resistance for each regenerator half is then lumped into a flow resistance co- 
efficient. The assumption is illustrated in Figure 6 which shows the three modes among 
which we wish to calculate the gas mass flow rate. The appropriate flow resistance coefficients 
and pressures are labeled. 


The mass flow rate into and out of either end of the regenerator is then calculated from 


/ 'p. - pj T, 

B V K Rj 


( for the ambient end 
of the hot or 
cold regenerator 


(39) 


Tp' 


Cl: 



- P.I P/R 
KR, 


for the hot or 
(cold end of the 
regenerators 


(40) 


where the constants KRj, KR^ are calculated from experimental data (e.g., as contained in Refer- 
ence 3) for the particular regenerator, and the sign in front of the radical is determined from an 
appropriate sign convention.* 


Now, in order to facilitate solving Equations (33) and (35) we make the following three approxi- 
mations in these two equations: 


*Details of the sign convention for Equations (39) and (40) c'e given in Reference 4. The signs in front of 
these equations ore of no significance in this discussion. 
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( 41 ) 


ni 
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X 


L 


RO 


M , 



(42) 


and 


(43) 


where again the signs in front of the terms are from a compatable sign convention. In effect, 
Equations (41) and (42) represent an average value for the regeneration mass flow rate and mass 
flow rate gradient, respectively. 

F. Internal Mass Flow Rates 

With the lumping of the regenerators’ flow resistance, the mass flow rates into and out of 
each cylinder can be calculated. Again, we assume the regenerator gases to be at an average 
pressure and temperature given by Equations (37) and (38), and we calculate the flow rate to 
a point at this temperature and pressure. Experimental data are used to obtain the flow resist- 
ance coefficient for each half of the regenerator. 

The internal mass flow rates can then be calculated by the equations 


and 




PCRO' I 


CRG 


^KR2 




PcRo' PcrgATcrg 

^R1 



2H 



(44) 

(45) 


(46) 

(47) 


G. Leakage and Friction 

The gas leakage mass flow rates around the cold and hot displacer are calculated from the 
laminar Couette flow solution 
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(48) 


M, » 


R R 


1)^ 

I 2 : 


(’’a- 


The frictional force exerted by the gas on each displacer can then be calculated via the 
laminar flow solution to the Navier Stokes equations; 




V 






(49) 


The heat transferred to or from the cylinder by the leakage gas is then approximated as 


C,. C^c - 2-a) (50) 

and 

Okl V<2p^T„-Ta). (51) 

where again we are assuming that the gas is heated or cooled to the gas temperatures in the 
cylinders to which it is flowing. 

H. Heat Transfer, Work, and Performance Computations 

The heat transfer rate to the hot cylinder is given by 

where the (T^^^^ - T^^^ ) term is the heat lost directly from the hot cylinder to the ambient 
cylinder (i.e., either by radiation or a direct conduction path). One-half of the heat transferred 
to or from the cylinder wall by the hot-to- ambient cylinder gas leakage is allocated to the hot 
cylinder. 

Similarly, the heat transfer rates to the cold and sump cylinder walls are given by 


and 


Oc ^CA "^2C (^^2C “ 


+ I _ 6 


CL 


(53) 
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^^AW “ "^A^ ’ ^CW ^ ^AW “ AW ‘ 


- %,, (-p < r„ - T„,, C„ ,T„ _ T^,) . ^6„, . 1 6^^, 


( 54 ) 


The average heat transfer rate over one revolution of the crankshaft is obtained by solving 
Equations (52) - (54) for the total heat transferred to each cylinder over one revolution, and 
dividing by the cycle period. Hence, 


0 


H 



(55) 


and 



(56) 


(57) 


The work that must be supplied to each displacer is the sum of the "pdV" work plus the work 
required to overcome the viscous drag of the gas between the displacers and cylinder walls. The 
work rate, then, is given by 


w K 


Vc( + 


\)C 


Pc) 


V- 


‘DH 




Pa V. 

A A 


(58) 


where and are computed from Equation (49) for the appropriate cylinder, with the sign de- 
termined from a compatable sign convention. 

The average work rate over one cycle is obtained by solving Equation (58) for the total work 
and dividing by the cycle period. Hence, 


WK 


WK 

T^CC 


(59) 


The coefficient of performance for the VM refrigerator is computed (by definition) from 


COP = 


^ 

(Q„ + :WKl) 


(60) 
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where the absolute value sign is required for the work rate term because WK is always negative as 
a result of the sign convention system employed. 

Another measure of performance used here is the percent of Carnot efficiency for a re- 
frigerator operating between T^, and T^. This is given by 


V. SYSTEM SOLUTION 


CAR 


COP 



( 61 ) 


A computer program (Reference 4) was developed which simultaneously solves the equations 
developed in sections IV A through IV H. This program employs a combination of finite difference 
and Runge Kutta numerical techniques to solve the system of equations as functions of crank angle 
. For most designs, the system consists of several hundred equations, primarily because the re- 
generators are subdivided axially into slabs (usually about 100 each); thus the energy and momen- 
tum equations must be solved for each slab at every . In addition, since the system represents a 
two-point boundary value problem, an iterative technique on must be employed. Typically, a 
step size of .004° is used with about 10 iterations and a run time of about 2-1/2 hours. 

Because of the large amount of computer time required, the program is not used for general 
VM refrigeration performance studies, but is primarily employed in the design of a specific 
Vuilleumier refrigerator. 

The input and output variables and constants for the program are listed in Appendix C. 


VI. RESULTS 

A. Data Correlation 


Very little design information or data on Vuilleumier refrigerators has been published. A 
design is presented, however, in Reference 1; a schematic of this machine, taken from Reference 
1, is shown in Figures 7A and 7B. The rather limited available data from testing this refrigerator 
are used to demonstrate program/test correlation. 

The results of the computer program runs for this refrigerator are shown in Table 2. The 
following points are noted with reference to this table and Figimes 8 through 10: 

(1) A difficulty was encountered in determining what effective cold regenerator mass to use 
for the program. The reason for this problem is that because of the low mass of the actual 
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Toble 2 


Program Results for Sample Refrigerator Design 


Run No. 

U 

SC 

5 

^CW 


Q 

H 

WK 

COP 

— 

CAR 

1 

.06 

-320" F 

1.67 

140.5 

5.75 

.0114 

.0326 

2 

.04 

-320"F 

1.11 

137.7 

5.7 

.0076 

.022 

3 

.0186 

-320° F 

- .03 


5.55 

- 

- 

4 

.04 

-357°F 

-2.8 

147.7 

7.65 

- 

- 

5 

.04 

-345°F 

-1.48 

142.1 

6.87 

- 


6 

.06 

-330° F 

.584 

141.8 

6.19 

.00394 

.0124 



Figure 8. Cooling Load vs. Effective Cold 
Regenerator Mass for Sample Refrigerator. 
Tc = -320^F. 


regen»,rator screening material, the obtainable cooling for a given cold cylinder temperature is 
sensitive to any additional regeneration which may occur at the inlet or outlet of the cold regener- 
ators. The various possible sites for regeneration at the inlet and outlet to the regenerator were 
thus investigated. It was determined that the stainless steel gas port insert at the base of the cold 
displacer (Figure 7) provides enough effective mass and heat transfer capability for regeneration. 
This effective mass was calculated by computing the effective depth of penetration of a sinusoidal 
tempera'.ure variation in the inside and outside surfaces of the insert (Reference 5). 
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Figure 9. Hot Cylinder Power Computations 
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Figure 10, Cooling Load Computations 
for Sample Refrigerator 


(2) Figure 8 shows the calculated cooling load vs, the effective cold regenerator mass for 
-320"* F. The lowest cooling load point includes only the screens as the effective mass, while 
the high point includes the effective mass of the whole inside and outside surfaces of the stainless 
insert. The middle point also includes as the effective mass the inside surface of the stainless in- 
sert, but only half the length of the outside surface (i.e., from the gas ports to just beyond the 
carbon rings). It was felt that because of the influence of the outside wall which is at the ambient 
wall temperature, the most realistic estimate of the effective regenerator n*ass is .04 lbs (i.e., 
the middle point in Figure 8.) Therefore, the correlations between theory and data which follow 
use this value for the cold regenerator mass. Finally, it is of interest to note that the convex shape 
of the curve in Figure 8 is expected, since the greater the cold regenerator mass, the less effect a 
change in the mass has on the obtainable cooling load. 


(3) Comparisons between the calculations and data are shown in Figures 9 and 10. 


Figure 9 shows that there is very good agreement between the data and the present analysis, 
especially when compared to the ideal analysis. Furthermore, most of the difference between the 
hot cylinder heat transfer rate of the data and of the analysis can be attributed to insulation and 
lead wire losses which are not included in the analysis. 

Figure 10 shows only two data points, with a linear interpolation between them. It is seen that 
the VM analysis results agree much more closely with the data than does the ideal model, especially 
at thiB higher cold-cylinder temperatures. For example, at Tc = 140°R (= 77K), the program pre- 
dicts an obtainable cooling load of 1.1 watts, as compared to 2.5 watts measured, while the ideal 
model predicts 12.8 watts. At the colder cylinder temperatures the agreement between the present 
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Figure 11. Pc “ Diagram for Sample Refrigerator 



Figure 12. Ph “V^ Diagram for Sample Refrigerator 


model and data is not as good, although still much 
closer than the ideal model. (Note chat a nega- 
tive cooling load means that at the correspond- 
ing cold cylinder temperature, the net heat would 

be transferred from the cold cylinder to the load.) The reason for the wider discrepancy at the 
lower cold cylinder temperatures might be the assumption of a constant matrix specific heat along 
the length of the cold regenerator. Obviously this assumption becomes less accurate as the differ- 
ence between the temperatures at the ends of the cold regenerator increases. 


(4) It is concluded that iood agreement exists between the program and the available data. 
Furthermore, it is felt that the agreement with respect to cooling load would be even better for a re- 
frigerator with a cold regenerator that was not quite as marginally designed (i.e., less sensitive to 
changes in regenerator mass). Finally, several suggestions are proposed in Section VII which should 
further improve the correlation for the cold end of the refrigerator. 

B. Illustrative Results of Program Utility 


The computed P-V diagrams for the hot and cold cylinders of the sample refrigerator are 
shown in Figures 11 and 12. The area enclosed by the cold-end curve divided by the period of 
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revolution represents the maximum possible cooling load for the refrigerator under the specified 
operating conditions, while the area enclosed by the hot side P-V diagram divided by the period of 
revolution represents the minimum possible power requirement for the refrigerator design. Hence, 
the feasibility of the design with respect to geometric considerations (i.e., hot cylinder size, cold 
cylinder size, dead volume, etc.) can be checked. 

The differences between the final calculated cooling level and hot cylinder power input com- 
puted from the program, and those calculated from the P-V diagrams are due regenerator inef- 
ficiencies, non-isothermal heat transfer from the gas to the cylinder walls, gas leakage around 
the displacers, and conduction losses. As an example, for the sample design 6^, and 6„ are 9.3 
watts and 70.5 watts respectively from the P-V diagrams, as compared to 1.11 watts and 137.7 watts 
for the final calculated values. It is obvious, therefore, that there are substantial losses from the 
aforementioned processes. 


Regenerator performance can be evaluated with the aid of Figures 13 and 14, which were 
composed from program data. Figure 13 shows the gas mass flc^v rate at the cold end of the cold 
regenerator as a function of crank position. The positive values are the mass flow rates of gas 



Figure 13. Mass Flow Rate of Gos at Cold End of Cold 
Regenerator vs, « for Sample RegenproXor. 



Figure 14. Motrix and Gas Temperatures of Cold 
of Cold Regenerator vs for Sample Refrigerate 
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leaving the regenerator and entering the cold cylinder. The negative values are the mass flow 
rates of gas leaving the cold cylinder and entering the cold legenerator. Figure 14 shows the 
matrix and gas temperatures at the extreme cold end of the cold regenerator as functions of 
crank angle. 

During steady state operation, as the crank turns from the zero degree position the helium 
flows through the regenerator into the cold cylinder; hence, the mass flow rate is positive in 
Figure 13. The gas, as it flows through the cold regenerator, is cooled to the temperatures shown 
in Figure 14. At the same time the temperature of the matrix at the cold end of the cold regener- 
ator increases as it is heated by the exiting gas. This process results in increasing the gas exit 
and matrix temperatures until the flow direction reverses at - 165° (Note the minor irregularity 
at the peak of the curves in Figure 14, in which the gas temperature curve crosses over the matrix 
temperature curve, while the gas mass flow rate is still positive. This irregularity appears to be 
due to a minor numerical instability at the area of the peak.) As the gas mass flow rate reverses, 
the gas temperature is identical to the cold cylinder gas temperature. The flow of this cold gas 
over the regenerator's cold end lowers the matrix temperature as shown in the figure. Following 
this, the cold cylinder gas temperature first increases and then decreases, with a corresponding 
matrix temperature as shown. Finally, as the flow rate direction reverses again at - 0°, this 
whole process starts over again. 

The regenerator's heat loss per cycle can be expressed by 


I 


^I.OSS I *'ic ^ 




T,c d . 


(62) 


(63) 


where Tj^ is the gas temperature shown in Figure 14 and is the mass flow rate shown in 
Figure 13. By means of Figures 13 and 14, therefore, the regenerator inefficiency can easily be 
calculated. Note, however, that this loss is coupled with the fluctuation of the cold cylinder gas 
temperature, since this fluctuation affects the temperature ^ving of the matrix. In this way, there- 
fore, the program can be used to design, or to improve existing designs, of regenerators. In the present 
case, it is obvious (without going through the above computation) from the large matrix temperature 
swing that an increase in the mass of the cold regenerator would significantly increase the obtain- 
able cooling loads. This conclusion was also reached in Section VI A. 

Another regenerator design variable which can be studied with the program is the gas pres- 
sure drop across the matrix. This is particularly important in designing the seals for a refrigerator 


22 



for long life-time operation. As an example, 
Figure 15 shows the cold and ambient cylinder 
pressures vs. crank angle for the sample re- 
frigerator, where the difference between the 
curves is the differentici pressure across me 
cold displacer seals. Similar curves for the hot 
cylinder show a maximum pressure drop on the 
order of 1 psi across the hot regenerator. 

Cylinder-to-gas heat transfer coefficients, 
gas leakage effects, and conduction losses can 
best be evaluated by making several computer 
runs for a specific refrigerator design, with 
several different values for the constants which 
represent these effects. 

VII. POSSIBLE PROGRAM IMPROVEMENTS 

The VM analysis program could be im- 
proved in the following ways: 

(1) The regenerator matrix specific heat 
could be allowed to vary with temperature ac- 
cording to input values or functions. 

(2) Regenerators composed of sections of 
simple program modifications. 



Figure 15. Ambient and Cold Cylinder Pressures 
vs. ^ for Sample Refrigerator 


different material could be handled by reasonably 
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Appendix A 
List of Symbols 

Ay hot or cold displacer cross-sectional area 
Ayj, cold displacer cross-sectional area 
Ay„ hot displacer cross-sectional area 
regenerator free flow area 
cold regenerator free flow area 
Ap„P hot regenerator free flow area 
b gap between displacer and cylinder 
CAR % of Carnot COP 

cold regenerator flow resistance coefficient defined by Equation (45) 
regenerator flow resistance coefficient defined by Equation (43) 

Cp gas specific heat at constant pressure 
Cpp regenerator matrix specific heat 
COP coefficient of performance 

CTAic initial gas temperature at ambient end of cold regenerator 
CTCIC initial gas temperature at cold end of cold regenerator 
Cy gas specific heat at constant volume 
D constant defined by Equation (36) 
d m infinitesimal element of gas in regenerator 
d infinitesimal element of gas entering or leaving hot cylinder 

d infinitesimal element of gas leakii^ around hot displacer between hot cylinder and am 

bient cylinder 
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r. constant defined l)y Equation (36) 

F hot or cold displacer frictional force 

F^ cold displacer frictional force 
hot displacer frictional force 
h regenerator gas enthalpy 

ambient cylinder- to~gas convective heat transfer coefficient times heat transfer area 
ambient heat exchanger H convective heat transfer coefficient times heat transfer area 
*^AHF .2 ambient heat e cchanger #2 convective heat transfer coefficient times heat transfer area 
cold cylinder -to-gas convective heat transfer coefficient times heat transfer area 
cold heat exchanger convective heat transfer coefficient times heat transfer area 
enthalpy of gas entering or leaving cold end of cold regenerator 
hot cylinder -to-gas convective heat transfer coefficient times heat transfer area 
hot heat exchanger convective heat transfer coefficient times heat transfer area 
H hot regenerator flow resistance coefficient defined by Equation (46) 

R 1 

H „ hot regenerator flow resistance coefficient defined by Equation (47) 

KR 2 

hj enthalpy of gas entering or leaving hot cylinder 

ii* convective heat transfer coefficient between regenerator matrix and gas 
li“ defined in Table 1 

HTAir initial gas temperature for ambient end of hot regenerator 
HTHic initial gas temperature for hot end of hot regenerator 

hj enthalpy of gas leakage entering or leaving hot cylinder via the hot displacer (Equation 
(16) 

I REG program print command 
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I RUN 
IRUNX 
ITT 
T 


K 


CA 


K 


cw 


K 


MA 


K„w 

KR, 

KRj 


L 


L 


0 


CD 


L 


HD 




M 

M 


tfi 


P. P 

B Mj/3 




program print command 
program print command 
number of iterations for program 
constant, equal to 778,16 ft-lb/BTU 

overall thermal conductance between cold cylinder gas and ambient cylinder gas, includ- 
ing conduction across the cold displacer 

thermal conductance from cold cylinder wall nirectly to ambient cylinder wall 

thermal conductance between hot cylinder gas and ambient cylinder gas, including con- 
duction across the hot displacer 

thermal conductance from hot cylinder wall directly to ambient cylinder wall 

regenerator flow resistance coefficient defined in Figure 6 

regenerator flow resistance coefficient defined in Figure 6 

cold cylinder connecting rod length 

regenerator length 

length of hot or cold displacer 

length of cold displacer 

length of hot displacer 

hot cylinder connecting rod length 

mass 

gas mass flow rate 
regenerator gas mass flow rate 
pressure 

regenerator mass flow rate defined by Equation (39) 
regenerator mass flow rate defined by Equation (40) 
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P .iverage hot or cold regenerator gas pressure defined by Equation (38) 
P‘ regenerator wetted perimeter 
hot or cold cylinder gas pressure 
o heat 

0 heat transfer rate 
r cranlt radius 
R gas constant 
Pp displacer radius 
cold displacer radius 
hot displacer radius 

HD 

RPM crank revolutions/minute 
SMr number of cold regenerator slabs 
SMU number of hot regenerator slabs 
t time 

T temperature 

hot or cold cylinder gas temperature 
T defined in Table 1 

HIN 

defined in Table 1 

hout 

u gas internal energy 
V volume 

w (defined in Table 1) 

WK work rate 

X coordinate along regenerator length 
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XDEG progrant print command 
XVAR program convergence constant 
WAR program convergence constant 

Greek Letter Symbols 
ri angle between hot and cold cylinder axes 
A crank angle step size 
A Cp/Cy for gas 
„ rotational speed (rad/sec) 

*1' hot and cold crank displacement angle (refer to Figure 5) 
cold-to- ambient leakage average gas viscosity 
hot -to- ambient leakage average gas viscosity 
hot- or cold-to- ambient leakage average gas viscosity 
^ gas density 

crank angle (zero position is top dead center of cold cylinder) 

Subscripts 

A ambient cylinder gas 

Ai gas entering or leaving ambient cylinder through heat exchanger #1 
A2 gas entering or leaving ambient section through heat exchanger #2 
AD ambient cylinder dead volume 
AW ambient wall 
c cold cylinder gas 
CD cold cylinder dead volume 
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CL leakage gas around cold cylinder 
CMAXO cold displaced 

CR cold regenerator dead volume 
CRG cold regenerator 
cw cold wall 
V. gas 

H hot cylinder gas 
HD hot cylinder dead volume 
HL leakage gas flow around hot displacer 
HMAXD hot displaced 

HR hot regenerator dead volume 
HRG hot regenerator 
HW hot wall 
L leakage 

or gas into or from ambient end of cold regenerator 
OH gas into or from ambient end of hot regererator 
1C gas flow from or into cold end of cold regenerator 
IH gas flow from or into hot end of hot regenerator 
R regenerator matrix 
RC cold regenerator matrix 

RG applies to hot or cold regenerator dead volume or length 

RH hot regenerator matrix 

2C gas flow from or into cold cylinder 
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2M gas flow from or into hot cylinder 
0 total mass of gas in VM refrigerator 

Superscript 

- average 





Appendix B 

Derivation of Cylinder Geometric Equations (7) through (9) 



Referring to Figure B1 , 


x=AR+AC-DC 


(Bl) 


- i + T - (-1 COS <p + r cos 6) 


(B2) 


= "C (1 - cos 5) + r (1 - cos fl). 


(B3) 


Now 


and hence 


r _ I 

sin 0 sin 


s in 0 = 


= s 


(l)sin« 


(B4) 

(B5) 

(B6) 
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Combining (B3) and (B6) 


X = f J 1 


cos 


r -% 

s i n" ' I s i n 


' + r (1 - cos 


But 


and 


Hence 


X - 2 r 

mu* 


’SlAX '^MAXD 


V _ 1 /, 

V 2 r \^ ~ 

*MAXD \ 


j^sin”* sin 


— (1 - cos ( 9 ). 


Now, adding the dead volume to this equation, we have 


V '^CMAXD It ,, 

Vc = — .<- (1 _ cos 


sin"* ^jsin ffj 


+ (1 - cos (9) ^ 


which is Equ?.tion (16) 

Equation (8) can be derived in exactly the same way, but with the substitution 


f?' = a + 6> - V. 


Equation (9) follows from a volume balance on one whole refrigerator. 


(B7) 


(B8) 


(B9) 


(BIO) 
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Appendix C 

Program Input and Output 

1 . Input 

The program input is shown on the next two pages, which are based on material from Refer- 
ence 4. 

2. The program output consists of the following information: 

a. Cylinder pressures as functions of 

b. Cylinder temperatures as functions of . 

c. All internal gas mass flow rates as functions of e. 

d. Cylinder volumes and volumetric time rates of change as functions of 6 . 

e. Cold and hot regenerator gas and matrix temperatures as functions of e . 

f. Hot cylinder power requirement, cold cylinder cooling load, motor power requirement 
and ambient heat rejection for the input conditions. 
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DATA 6 





















